Two-fluid model of a Bose-Einstein condensate in the cavity optomechanical regime 
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We analyze an atomic Bose-Einstein condensate trapped in a high-Q optical cavity driven by a 
feeble optical field, a situation formally analogous to the central paradigm of cavity optomechanics 
[Brennecke et al. (Science, 322, 235 (2008)]. We account for two-body interactions via a two-fluid 
model that retains the intuitive appeal of the optomechanical two-mode description. The Bogoliubov 
excitation spectrum of this system comprises a gapped upper branch and a lower branch that can 
include an unstable excitation mode. 



Cavity optomechanics offers much promise for appli- 
cations ranging from the detection of weak forces and 
small displacements to fundamental studies of the tran- 
sition between the quantum and the classical world [J]. 
One major stepping stone toward these objectives is to 
optically cool one or a few modes of vibration of a mov- 
ing mirror to near its quantum-mechanical ground state. 
This goal is likely to be achieved in the near future, in 
situations ranging from micro- and nano-mechanical sys- 
tems @, H, 0, 0, H 0, 0| to larg e mirrors as used in gravi- 
tational wave detectors @, Efl El, El . 

In parallel to these developments, recent studies of ul- 
tracold atomic gases trapped in optical cavities have pro- 
vided an alternative path to the study of cavity optome- 
chanics [lH, 14 , [l5| . Here the excitation of a collective 
mode of the cold gas plays the role of the vibrational 
mode of the mirror. Optomechanical studies using cold 
gases are appealing for a number of reasons. Most impor- 
tantly perhaps, methods for preparing the atoms in their 
motional ground state are well developed, allowing one 
to prepare the effective oscillator mode in its quantum 
ground state. Also, due to the high cooperativity of the 
cold atomic gas, one reaches the strong-coupling limit of 
cavity quantum-electrodynamics (QED) for mean intra- 
cavity photon numbers on the order of unity, offering a 
promising avenue to combined studies of cavity QED and 
cavity optomechanics. 

The optomechanical behavior of ultracold ato mic g ases 
was demonstrated in experiments by Gupta et al. [131 ] and 



by Murch et al. [lj|, who showed that the cavity field 



coupling to a collective center-of-mass excitation of the 
atoms results in oscillatory displacement of the g 
well as by Brennecke et al. [l5[ , who studied the coupling 
between a density modulated Bosc-Einstcin condensate 
(BEC) and the cavity field. This group also developed 
a simple two-mode description of the BEC that clearly 
illustrates the formal analogy with the moving mirror 
system, and is the basis for an accurate description of 
the observed dynamics. 

This Letter extends that model to calculate the low 
energy excitation spectrum of a BEC coupled to a weak 
light field. One key element is to describe the condensate 
by a two- fluid model that accounts for mode-broadening 



due to two-body collisions, while retaining the intuitive 
appeal of the two-mode description. We calculate the 
two-branch excitation spectrum via a Bo goli ubov- type 
perturbative treatment of the collisions 3, 17|. The up- 
per branch of the spectrum is gapped, while the lower 
branch can become positive imaginary for small mo- 
menta, corresponding to an unstable excitation mode 
that grows exponentially in time. 

We consider a BEC trapped inside a Fabry-Perot cav- 
ity and interacting with a standing wave light-field gen- 
erated by a pump laser of wave number k parallel to the 
cavity axis. For simplicity, we restrict our calculation to 
one spatial dimension along the propagation direction of 
the pump laser, assume that the trap is soft enough so 
that its impact on the ground state of the BEC can be ig- 
nored, and work at T = 0. In this Letter we focus on the 
effects of both atom-light and atom-atom interactions on 
the static properties of the condensate. 

This system is described by the Hamiltonian H = 



and 



-Wap +^aa, where Hkm is the kinetic energy term, 



H ap = U Q a? a 



L/2 



L/2 



dx&(x) [cos 2 (fca;)] *(a 



(1) 



is the atom-field interaction, where L is the cavity length, 
H = 1, and Uo the single-photon light shift of the atoms. 
Finally, 7i aa describes two-body collisions in the familiar 
fashion, with U aa the atom-atom interaction strength. 

We proceed by expanding the Schrodinger field opera- 
tor of the condensate in terms of plane waves, 



(2) 



where 4> q (x) = -^=e tqx , b q is the corresponding boson 
annihilation operator, and the sum runs over all momenta 
q satisfying periodic boundary conditions. In general, the 
atom-field interaction results in the recoil of atoms by 
±2£k, where I is an integer. Typically many momentum 
side modes are excited, but for feeble fields we expect 
intuitively that only the ±2fc side modes are appreciably 
occupied. The field operator then reduces to 



= 4>o(x) b + 4>2k(x) b 2 k + 4>-2k{x) b-2k ■ 



(3) 
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It is this truncation that results in the analogy between 
the optically driven BEC and radiation pressure driven 
mirrors. 

The validity of the restricted plane wave basis is deter- 
mined by the ratio between the effective Rabi frequency 
Uqu c and the recoil energy E rcc = k 2 /2m, where m is the 
atomic mass. One can get a sense of its region of validity 
by numerically calculating the mean-held ground state of 
the BEC subject to the potential V(x) = U$n c cos 2 (fcx), 
where n c is the mean intracavity photon number, and 
evaluating the occupation probability Ph of states with 
£ > 1. When P\ w <C 1, the truncation ([3]) captures the 
qualitative features of the system. Using the parame- 
ters of Ref. [IH [IH , this calculation suggests that we are 
safely in that regime for n c < 25 (Pho — 0.08). Since 
the optical field couples the condensate to the q = ±2fc 
side modes symmetrically we expand with respect 

to symmetric ('cosine') and anti-symmetric ('sine') oper- 
ators, 

$(x) = 0o (x) b + y / 2/Lcos(2kx) c + i> 



'2/ Lsm(2kx) §o, 
(4) 

where c = ^= (b 2k + b-2k) an d s = ^= (b 2k - b_ 2 k) 
obey boson commutation relations. In the absence of 
optical fields and in the collisionless regime the many- 
atom ground state is a pure condensate where each atom 
occupies the bo mode. Switching on the feeble optical 
field couples the mode 6o to cq, but leaves §o unoccu- 
pied. This is the situation described by the two-mode 
(or single-mirror) optomechanical model of Ref. [l5| . [l^ 
Two-body collisions complicate the situation by scat- 
tering atoms from the modes q = {0,±2fc} into nearby 
momentum states. The momentum distribution then 
consists of three sub-distributions, each centered about 
one of q = {0,±2fc}. If these sub-distributions are nar- 
row enough that they do not appreciably overlap, then 
one may treat each of them as a distinct Bose gas. The 
symmetry of the atom-field interaction suggests again the 
introduction of 'cosine' and 'sine' operators 
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(5) 



where the condition \q\ < k ensures that we can treat 
the components of the BEC described by the modes 
|S g , Sg, c 9 |, as three distinct fluids. The widths of their 
momentum distributions depend only on the mean-field 
interaction energy gn, where g = (47rfi 2 a s ) /m is the two- 
body interaction strength. The other relevant energy 
scale is P roc , hence one can get a sense for the valid- 
ity of the three-fluid model from the Bogoliubov ground 
state occupation distribution as a function of the ratio 
r = gn/E YOC . For the conditions of Ref. 15| we have 
r ~ 0.67 and we find numerically that the three-fluid 
model is indeed valid. 



In the case where the intracavity field is in a Fock state 
|n c ), or for classical fields, we have that a — > yjn~ c , with 
n c the mean intracavity photon number. [20| Within the 
three-fluid description, the collisionless Hamiltonian be- 
comes, in dimensionless units where lengths are in units 
of k~ x and energies arc in units of P rcc , 



Ho = V2eN + 



q 2 b\b q 



{CqCq 



s\s q ) 



q 



i 



SqCq + CqSq'j 



Here N is the total particle number operator, 



(hU n c ) I (2V2X 



(0) 



(7) 



and the sums are understood to span —1 < q < 1. 

The first term in H. is an unimportant energy shift, the 
second term is the kinetic energy, the third term describes 
the optical coupling of the modes centered around q = 
(b q ) to the 'cosine'-modes, and the fourth term accounts 
for the kinetic coupling of the 'sine' and 'cosine' modes. 
The effect of collisions is implicit in the inclusion of q =/= 
terms. 

The reduction from a three-fluid to a two-fluid descrip- 
tion proceeds by neglecting the last term in the Hamil- 
tonian ([6]) . The validity of this step depends on the ratio 
between the couplings in equation ([6]). Specifically, for 
4<7/e <C 1, the true ground state has all qualitative fea- 
tures of the two-fluid ground state, the sine mode con- 
tributing only a small perturbation. The two-fluid ap- 
proximation is particularly good for low momenta and 
n c relatively large, but still sufficiently small that the 
population of higher momentum side modes remains neg- 
ligible. Its validity can also be readily judged graphically, 
by comparing the collisionless spectra of the two-fluid ap- 
proximation of H.o to that of equation ^ . 

We mentioned already that the Bogoliubov spectrum 
of the two-fluid system consists of two branches. We 
concentrate in the following mostly on the lower branch 
of the spectrum. Carrying out an analysis using the 
experimental parameters of Brcnnccke et al. [TBI ], and 
n c = 10 we found that the inclusion of the 'sine' mode 
results in a negligible energy difference for q < 0.1, while 
E-2F — P3F — 0.2P rcc at q = 0.5. A related test is to cal- 
culate the occupation probability of the sine-mode in the 
lower branch of the three-fluid model. This probability 
is virtually zero for q < 0.2, and is ~ 0.025 at q = 0.5. 

The two-fluid version of the Hamiltonian ([6]) is readily 
diagonalized in terms of linear combinations of b q and c q . 
Carrying out straightforward algebra we find 



Ho = V2eN + [e/9 («) + e 7 (?) l\% 



(8) 



where 



e/3 



(q) = 2- v / 4 + £ 2 +g 2 ; e 7 (q) = 2+y / A + e 2 +q 2 . (9) 
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Figure 1: Dotted line: gn/hUon c , as a function of n c . Solid 
line: Unstable region of the two- fluid system - the lower 
branch of the Bogoliubov spectrum is unstable for values of q 
below that line. In the region where a perturbative treatment 
of the collisions is valid, the size of the unstable region shrinks 
as n c increases. Parameters as in Ref. [l5| (see also endnote). 



The bosonic "dressed basis" operators are related to the 
bare operators by 



U —v 

V /I 



(10) 



2 



£ 2 , and e° = 



where pL = e°/^/e° 2 + e 2 , v = y/ 
e 7 (0). 

A Bogoliubov-type perturbative treatment of two- 
body collisions is appropriate provided that gn is small 
compared to MJon c . Figure [T] plots (gn) / (HUon c ) as a 
function of n c . and indicates that already for very mod- 
est values of n c the optical interaction dominates over the 
collisions. In this regime, we apply the Bogoliubov ap- 
proach to the ground dressed state, that is, fto — ► y/No, 
where No is the corresponding macroscopic population. 
Then we keep only terms of quadratic order or less in the 
operators j3 q and \, and follow the usual procedure (l6l | 
to arrive at a linear eigenvalue equation for the Bogoli- 
ubov excitation spectrum. 

Labeling the diagonal operator for this approximation 
of H by 

Tq = Uqflq + V q {3l q + W q % + Z q rf_ q , (11) 

we find that the corresponding eigenvalue equation is 



(M 



0. 



(12) 



where V g = [u q v q w q z q \ , I is the identity matrix, and 

(13) 



M = 



K x J ir 2J 2 r J 2 r 
-Jir —Ki -J 2 r -2J 2 r 
2J 2 r J 2 r K 2 J 3 r 



-J 2 r ~2,hr -,hr -K 2 



The matrix elements of M are K\ 
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(2J 3 -,h)r + 2^/T+tf + 
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Figure 2: Coupling parameters Ji, J2, and J3 as a function of 
n c . The labeled ticks on the right-hand side of the frame de- 
note the J-parameter values in the limit n c — > 00. Parameters 
as in Fig. 1. 



\ (3z/ 4 + 12;/ V + 2//) , J2 = \ (3^V - 4^ 3 ) , and 
J3 = I [2v A — 3v 2 fi 2 + 2//) , indicating that M is char- 
acterized by the three coupling constants {Ji, J 2 ,Jz}. 
J\gn and J 3 gn can be thought of as the mean-field inter- 
action energy per-particlc corresponding to two-body col- 
lisions of atoms initially occupying the lower and upper 
branch of the dressed spectrum, respectively, and J 2 gn 
sets the per-particle energy scale for the interaction be- 
tween the two fluids. For n c > 0, J 2 is always negative, 
indicative of an attractive interaction. 

Figure [2] shows these coupling constants as a function 
of n c . For n c — * one recovers the familiar Bogoliubov 
spectrum of a weakly interacting Bose gas. For n c — > 00, 
Ji = ^p, J 2 = — g and J3 = I and the lower branch of the 
excitation spectrum assumes the traditional Bogoliubov 
form characterized by the mean-field interaction energy 
per particle ^fgn — however, as already discussed, the 
two- fluid model breaks down in that limit. Still, Fig. [5] 
shows that the J-parameters approach their limiting val- 
ues even for photon numbers well within the optome- 
chanical regime, which demonstrates the relevance of the 
n c — > 00 limit. 

The upper branch of the spectrum corresponds to the 
dressed mode whose main contribution stems from the re- 
coiled component of the condensate (the 'cosine mode'), 
while the lower branch has a similar correspondence to 
the 'condensate mode'. For vanishingly small intracav- 
ity fields, the upper branch is separated from the lower 
branch by 4£' rec , as expected. More interesting is the 
lower branch displayed in Fig. [3l which is characterized 
by a region of small momenta where the spectrum be- 
comes positive imaginary, indicative of an instability that 
grows exponentially in time. This instability finds its ori- 
gin in the attractive interaction between the two fluids, 
and its existence means that the two-fluid ground state 
is not the true ground state. For larger values of q, the 
lower-branch spectrum resembles a conventional Bogoli- 
ubov spectrum, and as n c increases, varies continuously 
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Figure 3: Lower branch of the Bogoliubov spectrum of the 
two-fluid system. For small values of q it is positive imaginary, 
corresponding to an excitation mode that grows exponentially 
in time. The instability disappears for n c — > (n c — > oo), in 
which case we recover the conventional spectrum of a scalar 
condensate characterized by the mean interaction energy gn 
(17pn/8). Parameters as in Fig. 1, and E rcc = h 2 k 2 /2m is 
the recoil energy. 



between the limits n c = (corresponding to interaction 
strength gn), and n c — > oo (corresponding to l7gn/8). 

The size of the unstable region depends on the strength 
of the atom-field interaction relative to the two-body 
scattering. Figure Q] shows the size of the unstable do- 
main as a function of n c . Below the solid line the lower 
branch of the Bogoliubov excitation spectrum is purely 
imaginary with positive magnitude. Figure [1] shows that 
this region grows in the limit where perturbation theory 
ceases to be valid. For n c > 5, the unstable region shrinks 
as n r increases, and vanishes in the limit ttt^ > oo. 

One important observation is that the Taylor series 
with respect to q (> 0) for the lower branch of the spec- 
trum only includes terms of even parity. In contrast, the 
corresponding Taylor series for the conventional Bogoli- 
ubov spectrum has only odd terms. As is shown in Fig.[3j 
this distinction is most important for small q where the 
lowest order terms dominate, and the conventional Bo- 
goliubov spectrum is approximately linear. 

In summary we have investigated the low energy exci- 
tation spectrum of a BEC in the "optomechanical region" 
where there is a clear analogy between the dynamics of 
the collective excitations of the BEC, and the motion 
of a radiation driven moving mirror. We have extended 
the two-mode model that describes a collisionless BEC 
to a two-fluid model that accounts for two-body colli- 
sions, and have determined their effect perturbativcly. 
We found that the atom-atom interactions are dressed by 



the light-field - the most striking effect being the effec- 
tive attraction between the different fluids. This attrac- 
tion can lead to an unstable ground state characterized 
by exponentially growing excitations at small momenta. 
At larger momenta the low energy spectrum resembles a 
conventional Bogoliubov spectrum with enhanced atom- 
atom interactions. Further work will include an analysis 
of a "multiple mirror" situation, as well as a full quantum 
description of the optical field. 



DSG acknowledges stimulating conversations with Ju- 
lia Meyer. This work is supported in part by the US 
Office of Naval Research, by the National Science Foun- 
dation, and by the US Army Research Office. 



[2] 
[3] 

[1 

[5, 

[6 

[r 

[9 
[10 

[11 
[12 
[13 
[14 
[15 
[16 

fi- 
lls 



[19] 



[20] 



T. J. Kippenberg and K. J. Vahala, Opt. Express 15, 
17172 (2007). 

M. D. LaHaye et al, Science 304, 74 (2004). 

D. Kleckner and D. Bouwmeester, Nature (London) 444, 

75 (2006). 

C. H. Metzger andK. Karral, Nature (London) 432, 1002 
(2004). 

K. L. Ekinci, Y. T. Yang, and M. L. Roukes, Journal of 
Applied Physics 95, 2682 (2004). 

A. Schliesser et al Nat. Phys. 5, 509 (2009). 
S. Gigan et al. Nature 444, 67 (2006). 

J. D. Thompson et al, Nature 452, 72 (2008). 
T. Corbitt et al, Phys. Rev. Lett. 98, 150802 (2007). 
M. Anderlini, F. Marino, and F. Marin, Phys. Rev. D 80, 
013001 (2009). 

T. Corbitt et al, Phys. Rev. Lett. 99, 160801 (2007). 

B. Abbott et al, N. J. Phys. 11, 073032 (2009) 

S. Gupta et al, Phys. Rev. Lett. 99, 213601 (2007). 
K. Murch et al, Nat. Phys. 4, 561 (2008). 
F. Brennecke et al, Science 322, 235 (2008). 

C. J. Pethick and H. Smith, Bose-Emstein Condensation 
in Dilute Gases (Cambridge University Press, 2002). 

N. N. Bogoliubov, Journal of Physics (USSR) 11, 23 
(1947). 

Experimental parameters for 87 Rb are: two-body col- 
lision strength g ~ 5.6 x 10 -51 J • m 3 , atomic density 
n = 3 x 10 20 m~ 3 , laser wavelength Xl ~ 780 nm, single- 
photon light shift U = 27T x 3.7 kHz 
This model can easily be extended to situations where 
higher-order momentum modes are excited by photon 
recoil, generalizing Eq. (O to include these modes. The 
situation becomes then analogous to an optomechanical 
system with £ mirrors coupled by the intracavity optical 
field. 

For arbitrary quantized fields one can diagonalize the 
collisionless Hamiltonian in analogy with the familiar 
dressed states approach of quantum optics. 



